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We present an introductory lecture on general features of cosmic rays, 
for non-experts, and some recent developments concerning cascade 
equations for air shower developments. 
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1 Discovery 

The history of cosmic rays goes back to Coulomb, who realized that a suspended 
charged metal sphere looses charge. This fact has been unexplained for a long 
time, until CTR Wilson (1900) found that this is due to the ionization of the 
surrounding air. 

In order to learn about the origin of this ionization, Wulf (1910) climbed the 
Eifel tower and found that the ionization decreases with increasing altitude, so 
the ionization seemed to be linked to ground. 

To further investigate this question, Victor Hess (1912) made a balloon flight 
up to an altitude of 5 km. He realized that the ionization increases again above 
1 km, which means that the radiation is coming from above. These findings are 
generally considered to be the discovery of cosmic rays. 

Pierre Auger (1932) finally detected coincidences between counters up to 
300 m apart, he indentified in this way the first extended air showers (EAS). 

Today we know: When a high energy proton or nucleus hits an air molecule, 
it creates many hadrons, which make further collisions, creating more particles, 
and so on, leading to a so-called hadronic cascade. 

Certain hadrons (mainly n°) decay into photons, each of them initiating an 
electromagnetic cascade (of photons and electrons). 
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2 Characteristics of Cosmic Rays 



In fig. [H we show the spectrum of cosmic rays as measured over the years by 
many experiments. One observes a surprisingly simple structure, the spectrum 
is essentially a power law E ~ 7 , with an exponent 2.7 below 3 10 15 eV and 3 above. 
The position where the slope changes is usually referred to as the knee. 

In particular interesting is the high energy end of the spectrum: Protons 
faster than 10 20 eV interact with the cosmic microwave background (making A 
resonances) , so they cannot come from far away; one expects therefore a galactic 
source. On the other hand the Larmor radius suggests an extragalactic source... 
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Figure 1 : The spectrum of cosmic rays 



3 Sources 

What are the sources of the high energy cosmic rays? Two scenarios are consid- 
ered [1,2]: 

• Top-Down scenario, which amounts to the decay of some extremely heavy 
particle 

• Bottom-Up scenario, where particles are accelerated by conventional elec- 
tromagnetic forces. Electric fields are difficult to find, so probably one has 
to deal with magnetic fields (Fermi) 

Fermi proposed that particles may be accelerated by entering and leaving moving 
magnetic clouds (like tennis balls), which leads to a gain of energy of AE/E = 
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4/3 f3 2 (second order), or particles are accelerated by entering and leaving moving 
shock fronts, with AE/E = 4/3 (3 (first order). 

In any case, we have AE = aE. Considering n iterations, one finds 

E = E (l + a) n . 

Supposing an escape probability e, we find for the number of particles 

N = N (l-e) n . 

So we get indeed an exponential spectrum, 

N / E \ ln ( 1 ^ £ )/ ln ( 1+a ) 
iVo = \E~oJ 

Possible accelerators are 

• Supernova shock waves < 10 15 eV. 

• Jets in AGNs: the jets create shocks into the intergalactic medium. 

• Rotating neutron stars (ms spin rates, B = 10 13 G): Dynamo effect ^10 20 V. 

4 Air Showers 

The first step is a primary hadronic interactions of a proton or a nucleus with 
air producing many hadrons. Then these secondary hadrons interact with air, 
giving rise again to hadron production, and so on. This is the so-called hadronic 
cascade. Hadrons may decay before suffering an interaction, for example: 

7T° -> 77 (98.8%) 

7T ± -> / 2 ± + u(u) (100%) 

K ± -> fi ± + u(u) (63%). 

Most muons reach the ground (at high energies) . 

Each 7 initiates an electromagnetic cascade (photons and electrons). In an 
EM cascade, photons interact via 

• Compton scattering with electrons from an atom (low E), 

• Photoelectric effect: the photon is absorbed by an atom, its energy is used 
to liberate an electron (low E). 

• Pair production: the photon interacts with the Coulomb field of a nucleus 
to produce a electron-positron pair (high E). 

Electrons interact via 
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• Ionization energy loss. 

• Bremsstrahlung: the electron interacts with the Coulomb field of a nucleus 
to produce a slowed down electron and plus the radiated photon. 

At high energy, the most important processes are pair production and Bremsstrahlung: 

e — > e + 7 

7 — > e + e, 

so in any case the number of particles doubles for each iteration. 

The splitting continues till some critical energy E c is reached (after this the 
particles get absorbed or decay only) 

We first discuss a toy model (Heitler 1944): 

• At each step the energy is given in equal parts to the two daughters. 

• A branching happens always after a collision length A, so the number of 
branchings at depth X is X/\. 

The number of particles and the energy per particle are 

N(X) = 2 X/X , E(X) = E /N. 
The shower maximum is at E = E c , so 

N(X max ) = 2 X — / A = ^ 

So we find the important general results : 

N max oc E , 
X max oc In E . 

The standard technique to calculate air shower developments is the Monte 
Carlo technique [3]: Each particle may decay or interact, whatever happens ear- 
lier, at the corresponding depth; new particles from decay or interaction are 
written to stack and treated subsequently. 

This technique is easy to implement, but very slow at high energies, since 
N oc E and even impossible at 10 20 eV. 

An alternative is the numerical solution of so-called cascade equations (CE) 
[4-7], which is at high energies considerably faster than the MC approach. 

5 Characteristics of Air Showers 

In the following we show results from the CE calculations (=Conex) for proton 
induced showers (if not mentioned otherwise), and some comparisons with the 
CORSIKA Monte Carlo calculations [3] . The cascade equations will be discussed 
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Figure 2: Produced hadrons at 10 A 14 eV. 




Figure 3: Produced leptons, photons at 10 A 14 eV. 



in the following chapters. The final aim being the solution of the full threedimen- 
sional cascade equations, we show here only results concerning the X (depth) 
dependence of particle numbers, based on the ID version of the approach, which 
makes use of the small angle approximation. There is no information about the 
lateral distributions. 
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Figure 4: Produced hadrons at 10 A 16 eV 
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Figure 5: Produced leptons, photons at 10 A 16 eV 



6 Hadronic ID cascade equations 

Here we discuss the equations for the energy spectra h n (E, X) of hadrons of type 
n at depth X, where the depth is defined as 
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X = p(x)dx, 



with p being the density of air. The decrease of h n (E, X) due to collisions is 

dhri. h-n. 



dX \ n 

where the mean inelastic free path X n is 

^air 



A, 



(n) ' 



The decay rate in the particle cm. system is dh n /dr = —h n /r n , so 

dhn m n p~ l (X) 



dX ct„E 



h 
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Figure 6: Comparison Corsika-Conex at 10 A 18 eV. 




The equation for hadron density, considering collisions and decay (only loss 
terms) is 



dh n (E,X) 
dX 



-h n (E,X)(- + ^^l 

\ n CT n E 



The full equation, including gain from higher energies, and a source term (gen- 
eralizing simple initial conditions) is 



dh n {E,X) 
dX 



-K.(E,X) 



1 m nP -\X) 



IK(E) ' cr n E 
h m (E',X) 



W mn (E',E) m m p-\X) 



^m{E' 

+S n (E,X). 



CT m E' 
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Figure 8: Shower maximum . 




Figure 9: Individual showers at 10 A 16 eV . 



where W mn (E',E) and D mn (E',E) are inclusive spectra of secondaries of type n 
and energy E produced in interactions [W) or decays [D) of primaries of type m 
and energy E'. 

To obtain a system of differential equations, one proceeds to an energy dis- 
cretization as 



which gives 



E = E ■ ■ d 



dhnijX) 

OX 



-h ni (X) 



m n p 1 (X) 

CT„Ei 



~\ + 



,ji ™mp \X) 
d~ m Ej 



-S n i(X), 
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with h ni (X) =h n (E i ,X), \ ni = \ n {Ei), and 



= [ ' f(^r)W mn (E 3 ,E)dE+ [ g (JL)W mn (E v E)dE, 

JEi^ JEi 



D ™ = J e ' f(§;) D rnn(E 3 ,E)dE+ [ g(^)D mn (E 3 ,E)dE, 



with 

f(x) = ^±, g(x) = l-f(x), 
c — 1 

in order to assure energy and particle number conservation. 
The solution of the homogeneous equation is 



h ni (X) = h ni (X ) exp 



(1 - W: n ) (X - X Q ) m n (L(X) - L(X )) 



\<n.i CT n E' t 



with L' = p 1 . The solution of the full equation is 

(1 - Wt) (X - X Q ) m n (L(X) - L(X )) 



h ni (X) = h ni (X Q ) exp 



v f ' r i>' ; ;;„ , n « m mP -\x> 



+ / EEtf 

Jx y m j =i+ i 



+ D 31 



+ S ni (X') 



exp [ (I - W^ n ) (X - X') _ m n (L(X)-L(X>)) \ x , 



7 Electromagnetic cascade equations 

We take the longitudinal coordinate to be z, the lateral coordinates x and y. The 
direction of a particle at a given point (x, y, z) may be written as 

(sin 9 cos 
sin 9 sin 
cos 9 

We define 

a = sin 9 cos , 
j3 = sin # sin . 

In the following we work with the coordinates a, (3, x, y, z of a particle and its 
energy E. We consider differential spectra of electrons, positrons, and gammas 

dN n 

l n (a,(3,x,y,z,E) 



da dp dx dy dE ' 



with n e {e + , e ,7}, for a given depth z. The energy E refers to the kinetic 
energy of electrons and total energy of photons. We suppose a source term 
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S n = S n (a, (3,x,y, z, E). In the following we will supress the arguments of the 
fuctions l n , S n etc, unless explicitely needed. The equations of electron-photon 
cascades may be written in the following form 

— * 

UV l n — L n [l] + Lcou\,n[ln] + S n , 

with 

_ - n d d n d 
u V = cos(J— + a—+f3—. 

oz dx ay 

and 

L n [l] = V / UE') W mn (E', E) dE', 

J E 

with the differential spectra 

W e - e -(E', E) = W hrems (E', E'-E) + W Moellei (E', E) Q{E - E cut ) - (a e -(E) - ^(3) 8(E' 

W e+e+ (E', E) = W hlcms (E', E'-E) + WWbha(£', E) Q(E - E cut ) - (a e +(E) - -^(3) 5(E' 

W e+e -(E',E) = W Bhah UE',E)&(E-E cut ), 

W ie -(E',E) = W pair (E',E) + W Comp ton(E',E'-E), 

W je+ (E',E) = W pair (E',E), 

W e ^(E',E) = W hrcms (E',E), 

W e+1 {E',E) = W hlcms (E',E) + W annih (E',E), 

W^(E',E) = W Compton (E',E)-a J (E)S(E'-E), 



and 



a n {E) = V [ E W nm (E,E')dE'. 
™ Jo 



The ionization energy loss coefficient is obviously zero for photons, /3 7 = 0, the 
same as the Coulomb term, Lcoui, 7 = 0. Otherwise 



£coui, e [y = J drjdC {l e (a + rj, (3 + (,...) - l e (a, (3, ...)} 



da Co u\(vX) 



di]d( 

where l e stands for l e - or l e + . An expansion gives 

Ua + ,,/» + { ... .)-««, A- 0= £ 

i, 3 
i + j > 

So 

f)2i-\-2j / 



da 2i d(3 2 i ' 
hJ 
i+j>0 
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with 



1 fj&j V 2 X 2j 1 d*Co*{&) 



So we finally have 

8 8 8 8 2iJr2 i] 



i + j > o 

In the small angle approximation one may integrate over the variables a,{3,x,y 
and one obtaines for 



l n (z,E) = j l(a, (3,x,y, z, E)dad(3dxdy 

the equation 

8 ~ 

which can be solved after energy discretisation without major problems. 

8 Summary (Cascade Equations) 

CONEX ID works already quite well, it gives practically the same results as 
the Monte Carlo calculations with Corsika, although the time needed for the 
computations is enormously reduced. Concerning the 3D version, there is work 
in progress. 
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